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We investigate a new way of using the quantum fluctuations in the vacuum as initial conditions for subsequent 
classical field dynamics. We show that with this method the field properly thermalizes at later times, whereas 
the method used previously leads to unphysical results. 



1. Initial conditions for a quench 

Investigation of initial value problems for quan- 
tum fields using numerical simulation are diffi- 
cult. The dynamics cannot be computed exactly 
and one has to use e.g. the classical approxima- 
tion, the large N or Hartree approximation. It 
is then not obvious how to specify initial condi- 
tions properly. Here we focus on this problem 
and investigate a new way to supply initial con- 
ditions for classical dynamics. We shall apply the 
same strategy to a simulation with the recently 
proposed |l[] Hartree ensemble approximation as 
well. 

As a test case, we consider a quench in a sim- 
ple 1 + 1 dimensional ip 4 model discretized on a 
lattice, with hamiltonian, 

H =£[3*£-3P*^x + §A£ + |A^. (1) 

x 

We model a sudden quench by assuming that at 
t = the sign of in the potential flips from 
positive to negative. At t = the system is in 
the vacuum state corresponding to the single- well 
potential (/i 2 > 0). Classically this would imply 
that the field is at rest, but triggered by (vac- 
uum) quantum fluctuations, the field modes with 
momentum \p\ < \p,\ will grow exponentially in a 
"spinodal decomposition" of the field. It is not 
obvious how one should supply quantum fluctua- 
tions as initial conditions for the classical fields. 
Until now one has followed the prescription j| 
that the classical field at t = should reproduce 
the one- and two-point functions evaluated in the 
quantum theory. 
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This implies that initial values for the (com- 
plex) Fourier amplitudes ip p and tt p of the classi- 
cal field and its velocity n = ip, must be drawn 
from a gaussian ensemble with probability distri- 
bution (using uj 0p — (p 2 + ^ 2 |) 1/2 ), 

P(tp,w) (x cxp[-^(|7r p | 2 +uj% p \(p p \ 2 )/2uj 0p ], (2) 
p 

An apparent problem with this prescription is 
that the initial vacuum fluctuations contribute to 
the energy density and the equilibrium tempera- 
ture that would correspond to the amount of en- 
ergy put into the field diverges when the cut-off 
is removed, T oc 1/a. One can only hope that in 
the early stages of the field evolution, this reser- 
voir of spurious energy in the short wavelength 
modes does not significantly affect the low-energy 
dynamics. 

These problems are absent in the new approach 
we shall describe next. During the spinodal de- 
composition, there may be a time window t_ < 
t < <+ in which ((jj p \ip p \ 2 ) 3> \ while the backre- 
action can still be ignored, ^A^ 4 <C i|/i 2 |</3 2 . If 
this is so, one can solve the quantum initial value 
problem analytically in the approximation A = 
for t < to with t- < to < t + . At time to we switch 
to a classical simulation since occupation num- 
bers (of the low momentum modes) have grown 
sufficiently large. As matching conditions we im- 
pose that the one- and two-point functions with 
subtracted vacuum fluctuations, arc reproduced in 
the ensemble of classical initial conditions. 

Using the A = solution, we can compute 
the one- and two-point functions, (ip p ) = (tt p ) = 
0, {\<p p \ 2 ) and (|7Tp| 2 ). The two-point functions 
contain contributions from vacuum fluctuations 
which we would like to subtract. Even though 
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the quench model at early times is clearly far out 
of equilibrium, we nonetheless use the following 
equilibrium relations to define the (time depen- 
dent) particle energy lo p and number density n p , 

(\ip p \ 2 ) = (n p + i)/wp, <|tt p | 2 ) = (n p + \)w p . (3) 

In this way we identify the contribution from vac- 
uum fluctuations by the | in these relations. To 
switch to classical dynamics, we specify a gaus- 
sian ensemble of classical fields at t = to with 
(\fp l \ 2 ) = n p /u) p and (|7r^| 2 ) = n p co p . Even 
though the analytical calculation gives a tail n p oc 
1/p 4 for large momenta, shall make the further 
simplification that n p = for the \p\ > /i. 

For comparison we shall also show results, still 
using classical dynamics, but using the customary 
J^] initial ensemble that includes vacuum fluctua- 
tions, n p — > n p + \ . In this case one might argue 
that the coupling constants need to be renormal- 
izcd, and that somehow the effect of the vacuum 
fluctuations should be subtracted. As this is a 
classical theory which is far from equilibrium, it 
is not clear how to implement such a renormal- 
ization in a consistent way. 

Finally we use the Hartree ensemble approxi- 
mation to compute the dynamics. Now we must 
specify initial conditions for the mean fields and 
for the mode functions. In order to be as close to 
classical dynamics as possible, we choose an en- 
semble of initial conditions, in which the mean 
fields have the same distribution as the classi- 
cal fields above at time to- Vacuum fluctuations 
are represented by the mode functions, which we 
choose as plane waves. We refer to ref. for 
further details of the Hartree ensemble method. 

2. Thermalization 

In order to study the effect of the different 
initial conditions, we monitor the thermalization 
process. Hence we plot the time dependence of 
the particle number density which, for classical 
fields, is defined as in (^) but without subtract- 
ing the |, n p = (d^KM 2 )) 1 / 2 . In Figs. ^ 
we show log(l + which is linear in u> for a 

Bose-Einstein distribution. We switch from the 
analytical solution to the numerical method at 
torn = 0.6; the results do not change when this 
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Figure 1. Particle number density using classi- 
cal dynamics with initial conditions that include 
vacuum fluctuations. The model parameters are: 
A/to 2 = 1/12, Lm = 102.4, 1/am = 5. 

time is chosen differently, provided it stays in the 
range .2^,tom^S2 . 

In Fig. [I] we use classical dynamics with the cus- 
tomary initial conditions including high momen- 
tum vacuum fluctuations, as in ref. [| At very 
early time tm = 5 — 10, one recognizes the unsta- 
ble modes from the large values of n at a^O.5. 
Also the vacuum fluctuations are clearly visible, 
since they imply that n w i or log(l + l/?i) w 1.1. 

At later times modes with larger momenta 
thcrmalize with a classical distribution, indicated 
by the curved line n p = T/lo p . Eventually, at 
tTO^3000 all modes have thermalized and the cor- 
responding temperature has become unphysically 
large. The temperature we obtain from fitting the 
(low oj) part of the spectrum, is plotted in Fig. 
^. After an initial decrease one clearly sees the 
unphysical rising of the temperature for tm^75. 

Next we use the new initial conditions, which 
do not include vacuum fluctuations. As can be 
seen in Fig. |^, the low-momentum modes at early 
times behave similarly to the previous case. But 
now the number density is exponentially sup- 
pressed for large momenta. At lager times we 
see that the distribution gradually moves towards 
a classical shape, but this happens much more 
slowly than in Fig. [l]. 

The lack of spurious vacuum energy in the sim- 
ulation of Fig. [2] also has a pronounced effect on 
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Figure 2. T/ie same as Fig. [1| 6m£ using initial 
conditions without vacuum fluctuations. 
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Figure 3. The same as Figs. [J and ^ but using 
the Hartree ensemble approximation. 

the temperature we obtain by fitting n against 
T/uj for low lo. Now the initial temperature drop 
continues, as energy continues to dissipate from 
the low momentum modes towards higher mo- 
menta. 

Finally we use the Hartree ensemble method to 
compute the field dynamics. For increasing time, 
more and more modes thermalize and align to- 
wards a Bose-Einstein distribution. At the largest 
time that we simulated, tm = 3000 the particle 
numbers follow a BE distribution over a signifi- 
cant range of energies. 

The temperature, obtained from fitting the 
low-momentum part of the spectrum with a BE- 
distribution, is shown in Fig. |]. At very early 
times, <^50, the behavior appears to be the 
same as in the previous classical simulations. At 
later times the temperature decreases very slowly, 
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Figure 4. Time dependence of the temperatures 
obtained from the particle densities at low lu of 
Figs. [J (top), (bottom) and^ (middle). 

much more slowly than with the classical dynam- 
ics. This is understandable because the particle 
number densities at high momenta remain expo- 
nentially suppressed. 

3. Conclusion 

Using vacuum fluctuations as seeds for subse- 
quent classical dynamics has to be done with care. 
The straightforward way, choosing an initial clas- 
sical ensemble that reproduces the quantum two- 
point functions, leads to misleading results at late 
times: the temperature of low-energy particles 
that equilibrate keeps rising and will eventually 
reach an unphysical value T oc 1/a. We have fur- 
ther demonstrated that the new prescription for 
initial conditions gives quite reasonable results, 
both when used with classical and with Hartree 
ensemble dynamics. 
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